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Abstract 

We review the fabrication and key transport properties of graphene double layers, consisting of two graphene 
monolayers placed in close proximity, independently contacted, and separated by an ultra-thin dielectric. 
We outline a simple band structure model relating the layer densities to the applied gate and inter-layer 
biases, and show that calculations and experimental results are in excellent agreement both at zero and in 
high magnetic fields. Coulomb drag measurements, which probe the electron-electron scattering between 
the two layers reveal two distinct regime: (i) diffusive drag at elevated temperatures, and (ii) mesoscopic 
fluctuation-dominated drag at low temperatures. We discuss the Coulomb drag results within the framework 
of existing theories. 
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1. Introduction 

Closely spaced double layer electron systems pos- 
sess an additional, layer degree of freedom, which in 
certain conditions stabilizes ground states with no 
counterpart in the single layer case. Notable exam- 
ples include fractional quantum Hall states (QHS) at 
even denominator fillings, such as v = 1/2 [HE] and 
z/ = l/4[3j|l],ora peculiar QHS at total filling factor 
v = 1 (layer filling factor 1/2) 0. The v = 1 QHS in 
interacting double layers displays striking transport 
properties such as enhanced inter-layer tunneling [6 
and counterflow superfluidity [7J [8], and has been 
likened to a BCS exciton condensate [9]. Dipolar su- 
perfluidity has been posited to also occur at zero mag- 
netic field [10] in spatially separated, closely spaced 
two-dimensional electron and hole systems, thanks to 
the pairing of carriers in opposite layers. Although 
remarkable progress has been made in the realization 
of high mobility electron- hole bilayers [Til IE] > an un- 
ambiguous signature of electron-hole pairing remains 



to be experimentally observed. The common thread 
in these phenomena is the inter-layer Coulomb inter- 
action being comparable in strength to the intra-layer 
interaction, leading to many-particle ground states 
involving the carriers of both layers. 

The emergence of graphene p~3j [14] [15] as an elec- 
tronic material has opened fascinating avenues in 
the study of the electron physics in reduced dimen- 
sions. Thanks to its atomically thin vertical di- 
mension, graphene allows separate two-dimensional 
electron systems to be brought in close proximity, 
at separations otherwise not accessible in other het- 
erostructures, and tantalizing theoretical predictions 
are based on this property [16j [17]. In light of 
these observations, it is of interest to explore elec- 
tron physics in closely spaced graphene double lay- 
ers. Here we discuss the fabrication, and key electron 
transport properties in this system, namely individ- 
ual layer resistivity and Coulomb drag. We introduce 
a model to describe the layer density dependence 
on gate and inter-layer bias, and show that calcula- 
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tions agree well with experimental results in zero and 
high magnetic fields. Coulomb drag measurements 
reveal two distinct regimes: (i) diffusive drag at ele- 
vated temperatures, and (ii) mesoscopic fluctuations- 
dominated drag at low temperatures. While we focus 
here on graphene double layers separated by a thin 
metal-oxide dielectric, a system with which the au- 
thors are most familiar with [T8l [T9] , we also note 
recent progress in graphene double layers separated 
by hexagonal boron nitride [20j [21] . 

2. Realization of graphene double layers 

The fabrication of independently contacted 
graphene double layers starts with the mechanical 
exfoliation from natural graphite of the bottom 
graphene layer onto a 280 nm thick Si02 dielectric, 
thermally grown on a highly doped Si substrate. 
Electron beam (e-beam) lithography, metal (Ni or 
Cr-Au) contact deposition followed by lift-off, and O2 
plasma etching are used to define a Hall bar device. 
The AI2O3 inter-layer dielectric is then deposited 
by atomic layer deposition (ALD), and using an 2 
nm thick evaporated Al film to nucleate the ALD 
growth. The total inter-layer dielectric thickness for 
the samples used our study ranges from 4 nm to 9 
nm. To fabricate the graphene top layer, a second 
monolayer graphene is mechanically exfoliated on 
a Si02/Si substrate. After spin-coating poly(metyl 
metacrylate) (PMMA) on the top layer and curing, 
the underlying Si02 substrate is etched with NaOH, 
and the top layer along with the alignment markers 
is detached with the PMMA membrane. The PMMA 
membrane is then aligned with the bottom layer 
device, and a Hall bar is subsequently defined on the 
top layer, completing the graphene double layer. 

We focus here on data collected from two samples, 
labeled 1 and 2, both with a d = 7.5 nm thick AI2O3 
inter-layer dielectric, and with an inter-layer resis- 
tance larger than 1 GQ. The layer mobilities are 
~ 10,000 cm 2 /V-s for both samples. The layer re- 
sistivtities are measured using small signal, low fre- 
quency lock-in techniques as function of back-gate 
bias (Vbg), and inter-layer bias (Vtl) applied on 
the top layer. The bottom layer is maintained at the 
ground (0 V) potential during measurements. The 
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Figure 1: (a) Schematic representation of a graphene double 
layer, consisting of two graphene monolayers separated by a 
thin AI2O3 dielectric. The two layer densities can be indepen- 
dently controlled using the back-gate (Vbg) bias along with 
the inter-layer (Vtl) bias applied on the top layer. Lower 
right: optical micrograph of a graphene double-layer device. 
The red (blue) contour marks the bottom (top) layer. The 
scale bar is 5 /im. (b),(c) Band diagram of a graphene dou- 
ble layer under an applied back-gate [panel (b)] or inter-layer 
[panel (c)] bias. Adapted from Ref. [19]. 

data discussed here are collected using a pumped 3 He 
refrigerator with a base temperature T = 0.4 K. 

To understand the layer resistivity dependence on 
gate and inter-layer bias, it is instructive to examine a 
band structure model which relates the applied Vbg 
and Vtl biases to the top (ny) and bottom (n#) 
layer densities [Figs. l(b,c)]. The applied Vbg can 
be written as the sum of the electrostatic potential 
drop across the bottom Si02 dielectric and the Fermi 
energy of the bottom layer: 

eV B o = e 2 (n B + n T )/C S i0 2 + E F (n B ) (1) 

E F (n) represents the Fermi energy of graphene rela- 
tive to the charge neutrality (Dirac) point at a carrier 
density n; n and E F (n) are positive (negative) for 
electrons (holes). Csi0 2 is the Si02 dielectric capaci- 
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Figure 2: (a) Carrier density contour plot of the bottom (left 
panel) and top (right panel) layer, calculated as a function of 
Vbg an d Vtl- (b) Longitudinal resistivity measured in the 
bottom (left panel) and top (right panel) layer in sample 1 at 
T = 0.4 K. Adapted from Ref. [19]. 



tance per unit area. Similarly, an applied Vtl can be 
written as the sum of the electrostatic potential drop 
across the AI2O3 dielectric, and the Fermi energy of 
the two layers: 

eV TL = E F (n B ) - [E F (n T ) + e 2 n T /C Al2 o 3 ] (2) 

A positive Vtl applied on the top layer induces elec- 
trons (holes) in the bottom (top) layer, which ex- 
plains the negative sign for the right hand side terms 
in Eq. (2). While the above equations implicitly as- 
sume that the two graphene layers are charge neutral 
at Vbg = Vtl = V, a finite doping of the two lay- 
ers can be included in the above model using additive 
constants to the left hand side in Eqs. (1,2). Alterna- 
tively, Vbg and Vtl can be referenced with respect 
to the bias values which render both layers charge 
neutral, a convention which we adopt in discussing 
our experimental results [22] . 



3. Magnetotransport properties of individual 
layers 

Figure 2(a) data show contour plots of and 
riT calculated according to Eqs. (1,2), with the fol- 
lowing dependence of the Fermi energy on density 
Ef(u) = hvpy/TT • n, consistent with masseless par- 
ticles with a Fermi velocity vp- The bottom Si02 
dielectric capacitance is Csi0 2 = 12 nF-cm -2 , a 
value determined experimentally by probing the ca- 
pacitance of metal pads located in proximity of the 
graphene double layer and confirmed by Hall mea- 
surements. The inter-layer capacitance value used 
here is Cai 2 o 3 = 340 nF-cm -2 [23], and the Fermi 
velocity vp = 1.15 x 10 8 cm/s. Figure 2(a) data show 
that the charge neutrality point of the bottom layer 
follows an almost linear dependence on Vbg an d Vtl, 
while the density of the top layer is controlled primar- 
ily by the inter-layer bias, and to a much lesser extent 
by the back-gate bias. To better understand these 
observations, let us first neglect Ep{n) in Eqs. (1,2): 
riT is controlled by Vt l only, and n b depends linearly 
on Vbg and Vtl, with Csio 2 an d Cai 2 o 3 respectively 
as proportionality coefficients. Although this behav- 
ior resembles Fig. 2(a) data, noticeable departures 
can be observed as a result of the reduced density of 
states in graphene, and consequently non- negligible 
Ep(n). First, np does depend on Vbg, implying an 
incomplete screening of the back-gate induced elec- 
tric field by the bottom layer, or equivalent ly a finite 
density of the states in the bottom layer, hence finite 
Ef{tib)' Second, the tlb = line shows a noticeable 
non-linearity as a function of Vbg an d Vtl near zero 
gate bias, a consequence of the non-negligible Fermi 
energy of the top layer. 

Figure 2(b) data show contour plots of the longitu- 
dinal resistivities pB (left panel) and pT (right panel) 
measured in the bottom and top layers, respectively, 
at a temperature T = 0.4 K. A comparison with Fig. 
2(a) data reveals a very good agreement between the 
experimental data and calculations, validating the 
model of Eqs. (1,2). Interestingly, this model can be 
used to make a direct measurement of the Fermi en- 
ergy in graphene, by employing one of the two layers 
as a resistively detected Kelvin probe [19 . Indeed, 
setting riT — in Eq. (2) yields eVrL = Ef(ub)' 
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Figure 3: (a) Contour plots of (left panel) and tit (right panel) calculated as a function of Vbg an d Vtl at B = 8 T. (b) 
pB (left panel) and pT (right panel) measured as a function of Vbg an d Vtl at B = 8 T and T = 0.4 K in sample 1. The data 
show QHSs marked by vanishing resistivity. The QHS filling factors are marked in both panels. 



The latter relation implies that the inter-layer bias 
required to bring the top layer to the charge neutral- 
ity point is equal to the Fermi energy of the bottom 
layer expressed in units of eV. Consequently, track- 
ing the top layer charge neutrality point in Fig. 2(b) 
(right panel) yields the bottom layer Fermi energy as 
a function of Vbg- 

We next address the layer density dependence on 
Vbg and Vtl in a perpendicular magnetic field (B). 
In the presence of a 5-field the electrons occupy Lan- 
dau levels (LL) with energies En = ±v Fy / 2heB\N\, 
where N is an integer denoting the LL index. Includ- 
ing the spin and valley degrees of freedom each LL has 
a degeneracy of 4eB/h. Quantum Hall states (QHSs) 
develop at filling factors v = ±4(7V-hl/2). Neglecting 
LL disorder-induced broadening, the Fermi energy 
dependence on density therefore writes Ep = E^, 
where N = Int[nh/AeB] and Int is the nearest in- 
teger function. Figure 3(a) data show contour plots 
of riB (left panel) and tlt (right panel) calculated as 
a function of Vbg and Vtl at B = 8 T using Eqs. 



(1,2), along with the same capacitance values used 
in Fig. 2(a): C Si o 2 = 12 nF-cm" 2 , C A i 2 o 3 = 340 
nF-cm -2 . In addition, a finite broadening is assumed 
for each LL using a Lorentzian density of states dis- 
tribution with a width 70 = 14 meV for the N = 
LL and 7^ = 6.5 meV for |JV| > [19]. 

Figure 3(b) data show the measured pB (left panel) 
and pT (right panel) as a function of Vbg and Vtl in 
sample 1, at T = 0.4 K. The data display alternat- 
ing regions of vanishing resistivity corresponding to 
quantum Hall states, and high resistivity regions cor- 
responding to half filled LLs. A comparison with Fig. 
3(b) data reveals a very good agreement of the charge 
neutrality point dependence on Vbg and Vtl for both 
layers, effectively validating Eqs. (1,2) model in high 
magnetic fields. Similar to Fig. 2(b) data, the Vtl 
value when the top layer is charge neutral is equal to 
the Fermi energy of the bottom layer. Consequently, 
the top layer charge neutrality line of Fig. 2(b) traces 
the bottom layer Fermi energy as a function of Vbg, 
and shows the staircase dependence corresponding to 



4 



242 K 


166 K 

• * 




1 1 1 1 1 

81 K 



n T (10 1 7crri ) 



-2 2 

n T (10 11 /cm 2 ) 



■ 



(£1) 



Figure 4: Contour plot of pn measured in sample 1 as a function of ng and ut, and at different T- values as indicated in each 
panel. The diffusive Coulomb drag shown here increases with T, is positive (negative) for opposite (identical) carrier polarity, 
and has a similar dependence on layer density in all four quadrants. 



discrete LLs. 

4. Coulomb drag of massless fermions in 
graphene 

An important attribute of a double layer system 
is the ability to probe the Coulomb drag between 
the two layers. A charge current (I Drive) flown in 
one (drive) layer results in a net momentum transfer 
to the opposite (drag) layer, thanks to the Coulomb 
interaction between electrons in the two layers. If 
no current is allowed to flow in the drag layer, a 
voltage Vorag builds up in order to counter the mo- 
mentum transfer. The drag resistivity is defined as 
Pd = {V Drag I ^ Drive) x (W/L), where L and W are 
the length and width of the sample area over which 
the drag voltage is measured. In effect, the drag resis- 
tivity is proportional to the scattering rate between 
electrons located in opposite layers. The dependence 
of pd on ub-, tit, and the spacing d between the 
two layers can provide insight into the interaction be- 
tween the two layers, as well as the ground state in 
individual layers. 

Coulomb drag between closely spaced carrier sys- 
tems has been probed in a variety of GaAs/AlGaAs 
heterostructures which, depending on the type of 



sample used, allowed experimental access to electron- 
electron drag [24j [25], hole- hole drag [26] , and 
electron- hole drag [23 [TTJ [12]. Assuming that the 
ground state in each layer is a Fermi liquid, and that 
the layer spacing is sufficiently large such that the 
Fermi wave- vector kp ^> dr 1 , the drag resistivity de- 
pends on the layer densities, temperature, and spac- 
ing as p D oc T 2 /(n 3 J 2 • n^ 2 • d 4 ). 

In Fig. 4 we show contour plots of pD as a func- 
tion of riB: and tit, measured in sample 1 at different 
temperatures. The data were measured by flowing 
current in the bottom layer, and measuring the drag 
voltage in the top layer as a function of Vbg and 
Vtl: which were subsequently converted into layer 
densities using Eqs. (1,2). We note that because 
the drag voltage is orders of magnitude smaller than 
the longitudinal voltage drop in the drive layer, very 
small inter-layer leakage currents can introduce arti- 
facts in the experimental data; care was taken here to 
ensure the measured drag voltage is not affected by 
inter-layer leakage. Figure 4 data explore pD in the 
four quadrants of the (n#, ut) plane, and reveals sev- 
eral noteworthy findings. First, pD is negative (posi- 
tive) for same (opposite) carrier polarity. Second, pD 
increases with increasing T. Third, pD vanishes at 
large ub,t values, and in the vicinity of the ub,t = 
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Figure 5: (a) Contour plot of pp. vs. (VbGiVtl) measured at T = 4.2 K in sample 2, using the top (bottom) layer as the 
drag (drive) layer. The data shows large mesoscopic fluctuations which follow the constant density lines of the top layer, (b) 
Contour plots of pp. vs. {ubiTlt) (left panel) and <jp vs. {ubiTlt) (right panel) corresponding to panel (a) data. 



points. Lastly, the pp> vs. (np^np) data in different 
quadrants are very similar, except for the sign change 
depending on the carrier type. As we explain below 
all these observations are consistent with Coulomb 
drag of masseless fermions where the ground state in 
each layer is a Fermi liquid. 

Several theoretical studies examined the Coulomb 
drag in graphene to date [28j El E3 El EI E3 EI] . 
In the weak-coupling limit, defined as kp d^> 1, the 
drag resistance scales as [28, 32 , 33| [34] 



h 



PD OC 
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(3) 



In the strong-coupling limit, defined as kp • d <C 1, 
pp. is expected to have a weaker dependence on layer 
density, and to be independent of the layer spacing 
as [33] E2 

h T 2 



p D oc -— • 



1/2 1/2 



(4) 



A close examination of Fig. 4 data reveals a In- 
dependence which is stronger than pp. oc T, but 



weaker than pp> ex T 2 , in the T-range probed here. 
A best fit to pd vs. n#,r along the diagonals tlb — 
fir = n, and tlb — —tit = n yields po oc n~ a , with 
an exponent a = 1.7 little dependent on tempera- 
ture. The a- value is lower than 3, expected in the 
weak-coupling limit, but larger than 1 expected in 
the strong-coupling limit. This can be readily un- 
derstood since kp • d varies between 0.4 and 1 in the 
density range where pp> is not vanishingly small. At 
low densities, where kp-d <C 1 disorder-induced pud- 
dles in the two layers lead to a decrease of pD towards 
zero. The experimental pp. dependence on n#, rip, 
and T is in good agreement with the expected theo- 
retical dependence of diffusive Coulomb drag in the 
Fermi liquid regime, although a quantitative match- 
ing between experiment and theory remains to be 
established. 

As the temperature is reduced below 50 K, the drag 
resistivity reveals a dramatic departure from the dif- 
fusive Coulomb drag of Fig. 4. In Fig. 5(a) we show 
a contour plot of pp> measured as a function of Vbg 
and VpL: at T = 4.2 K in sample 2, and using the 
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top (bottom) layer as the drag (drive) layer. Un- 
like the diffusive Coulomb drag of Fig. 4, the data 
of Fig. 5(a) show a pattern of large pB fluctuations 
centered around a zero average. The mesoscopic fluc- 
tuations increase in amplitude with reducing T, and 
completely obscure the diffusive drag at low temper- 
atures [18] . The pD mesoscopic fluctuations emerge 
as a result of phase coherent transport at low tem- 
perature, and represent the counterpart of universal 
conductance fluctuations in Coulomb drag [35 . Sim- 
ilar mesoscopic fluctuations have been observed in 
Coulomb drag probed between two-dimensional elec- 
tron systems in GaAs/AlGaAs heterostructures, al- 
beit at lower temperatures [36] . 

A comparison of Fig. 5(a) data with the carrier 
density contour plot of Fig. 2(a) calculations re- 
veal an interesting observation. The locus of con- 
stant pD in Fig. 5(a) corresponds to the lines of 
constant tlt in the (Vbg,Vtl) plane of Fig. 2(a). 
Alternatively stated, the pjj mesoscopic fluctuations 
track the drag layer constant density lines in the 
(Vbg,Vtl) plane. To better illustrate this observa- 
tion in Fig. 5(b) (left panel) we show a contour plot of 
pD vs. (n£,?7/r); these data represent Fig. 5(a) pB 
data with (Vbg,Vtl) axes converted into (n^,^) 
using Eqs. (1,2). The right panel of Fig. 5(b) (right 
panel) shows the drag conductivity <jb = Pd/(pb'Pt) 
vs. (n#,nT). We note that the gb fluctuation am- 
plitude has the same order of magnitude as e 2 /h. 
Figure 5(b) data manifestly show that the Coulomb 
drag mesoscopic fluctuations depend mainly on the 
drag layer density (nr), and are largely insensitive 
to the drive layer density (n#). This observation is 
manifestly at variance with the Onsager reciprocity 
relation, according to which interchanging the drag 
and drive layers should not affect the measured pB- 

Examination of Fig. 5(b) data reveals a regular, 
almost periodic pB vs. tit pattern. To better illus- 
trate this pattern, in Fig. 6 we show pB vs. tlt mea- 
sured at a fixed drive layer density, ub = —5 x 10 11 
cm -2 , and at T = 4.2 K. The data reveals an al- 
most periodic dependence of pB vs. tit reminiscent 
of Fabry- Perot interference [37]. A simple analysis 
can relate the top layer Fermi energy (AE F ^ T ) change 
corresponding to pB maxima to a cavity length (L), 
via AE Fj t = h • vf / L. A typical value for AEf,t 
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Figure 6: p]j vs. measured at a fixed bottom layer density 
riB = —5 x 10 11 cm -2 in sample 2, at T = 4.2 K and using 
the top (bottom) layer as the drag (drive) layer. The top axis 
represents the Fermi energy of the drag layer. The data show 
an almost periodic fluctuation pattern. 

of 10 meV deduced from Fig. 6 data, corresponds 
to L = 0.5 /im, a value comparable to the phase co- 
herence length in graphene on Si02 [38] at T = 4.2 
K. 

5. Conclusions 

In summary, we present a magnetotransport and 
Coulomb drag study in graphene double layers, con- 
sisting of two independently contacted graphene 
monolayers separated by a thin AI2O3 dielectric. The 
Coulomb drag probed in this system reveals two 
regimes: (i) diffusive drag at elevated (T > 50 K) 
temperatures, and (ii) mesoscopic fluctuations dom- 
inated drag at low temperatures. The temperature 
dependence of the diffusive drag is consistent with 
the Fermi liquid theory, while the density dependence 
suggests the layers are close to the strong-coupling 
regime [33, 34]. The Coulomb drag mesoscopic fluc- 
tuations observed at low temperature depend mainly 
on the drag layer density, and are largely insensitive 
to the drive layer density, an observation which is at 
variance with the Onsager reciprocity relation. 
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